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ABSTRACT 

We report final results from our 2.5 year infrared parallax program carried out with the European 
Southern Observatory 3.5m New Technology Telescope and the SOFI infrared camera. Our program 
targeted precision astrometric observations of ten T type brown dwarfs in the J band. Full astromctric 
solutions (including trigonometric parallaxes) for nine T dwarfs are provided along with proper motion 
solutions for a further object. We find that HgCdTe-based infrared cameras are capable of delivering 
precision differential astrometry. For T dwarfs, infrared observations are to be greatly preferred over the 
optical, both because they are so much brighter in the infrared, and because their prominent methane 
absorptions lead to similar effective wavelengths through the J-filter for both target and reference stars, 
which in turn results in a dramatic reduction in differential colour refraction effects. We describe a 
technique for robust bias estimation and linearity correction with the SOFI camera, along with an 
upper limit to the astrometric distortion of the SOFI optical train. Colour-magnitude and spectral- 
type-magnitude diagrams for both L and T dwarfs are presented which show complex and significant 
structure, with major import for luminosity function and mass function work on T dwarfs. Based on the 
width of the early L dwarf and late T dwarf colour magnitude diagrams, we conclude the brightening of 
early T dwarfs in the J passband (the "early T hump") is not an age effect, but due to the complexity 
of brown dwarf cooling curves. Finally, empirical estimates of the "turn on" magnitudes for methane 
absorption in field T dwarfs and in young stars clusters are provided. These make the interpretation of 
the T6 dwarf a Ori J053810. 1-023626 as a cOri member problematic. 

Subject headings: parallaxes - infrared: stars - infrared: brown dwarfs 



indicate that ~10-20% of objects observed in sufficient de- 
tail are found to be binary.) Luminosity functions based on 
currently available colour-magnitude relations will there- 
fore be problematic at best. Trigonometric parallaxes are 
therefore required to determine the actual luminosities of 
these objects and indicate whether they are single or bi- 
nary, so that more meaningful luminosity functions for T 
dwarfs can be constructed. 

2. PARALLAXES AND THE INFRARED 

Traditional parallax techniques based on photography 
are completely unable to target objects as faint and red 
as T dwarfs. CCD parallax work in the optical at the 
USNO, ESO and Palomar (Monet et al. 1992; Dahn et al. 
2002; Tinney 1996; Tinney et al. 1995; Tinney 1993) have 
shown that parallaxes can be obtained for objects as faint 
as 1=18-19 at distances <70pc. However, this still leaves 
the T dwarf class of objects (with I>21) unobservable. To 
date only a few of the very brightest and closest T dwarfs 
have proved tractable for CCD parallax work (Dahn et al. 
2002). 

Over the last two years, therefore, we have been extend- 
ing optical CCD astrometric techniques into the infrared, 
where the J<16 magnitudes of most of the detected T 
dwarfs make significant progress possible. Indeed, there 
are several reasons to prefer the infrared for high preci- 
sion astrometry. First, the effects of differential colour 

1 Based on observations obtained at the European Southern Observatory, Chile. Programmes 65.L-0061, 66.C-0404, 67.C-0029, 68.C-0004 & 
69.C-0044. 

2 Anglo-Australian Observatory, PO Box 296. Epping. 1710 Australia. cgtOaaoepp.aao.gov.au 

3 UC Los Angeles, 8371 Mathematical Sciences, CA. 90095. USA. adamaastro.ucla.edu 
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1. INTRODUCTION - METHANE T-TYPE BROWN DWARFS 

Numerous examples of the field counterparts to the ex- 
tremely cool methane brown dwarf Gl 229B (Nakajima et 
al. 1995) are now known (Strauss et al. 1999; Burgasser 
et al. 1999, 2000a,b; Leggett et al. 2000; Tsvetanov et al. 
2000; Cuby et al. 2000). These objects are now uniformly 
classified as "T dwarfs" (Burgasser et al. 2002b; Geballe 
et al. 2002) , and have such low photospheric temperatures 
(800-1300K), that their photospheres are dominated by 
the effects of dust and methane formation (Allard et al. 
2001), neither of which are amenable to simple modeling. 
The discovery of sizable numbers of T dwarfs, means that 
we are now in a position to use direct trigonometric par- 
allax observations to empirically determine the loci of T 
dwarf cooling curves, rather than relying on models. The 
discovery of several T dwarfs by SDSS with spectra bridg- 
ing the L and T spectral types (e.g. Leggett et al. 2002; 
Geballe et al. 2002) means we are also in a position to 
empirically determine where on these brown dwarf cooling 
curves the L-T transition occurs. 

Trigonometric parallaxes are also essential to under- 
standing the space density of T dwarfs. Luminosity func- 
tion estimates for T dwarfs (eg. Burgasser 2002a) are cur- 
rently based on limited parallaxes and assumptions about 
object binarity. (Recent programs targeting more L and T 
dwarfs (Martm, Brander & Basri 1999; Koerner et al. 1999; 
Reid et al. 2001; Burgasser et al. 2003a; Close et al. 2003), 
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refraction (the different amount of refraction the atmo- 
sphere produces in red target stars, compared to blue ref- 
erence stars, see Monet ct al. 1992) are reduced by working 
at longer wavelengths. Second, because T dwarfs suffer 
methane absorption at the red end of their J- and H-band 
spectra, their effective wavelengths through a J- or H-filtcr 
are much closer to that of a typical background reference 
star, than is the case in the optical. These effects combined 
mean that the stringent requirements on maintaining con- 
trol of observations at constant hour angles (at least for 
T dwarfs) is not present in the infrared (cf. Section 6.1). 
This considerably increases the flexibility and efficiency of 
infrared parallax observing, over the optical. Third, see- 
ing improves in the infrared, leading to smaller images, 
smaller amounts of differential seeing, and so higher astro- 
metric precision. And finally, T dwarfs show much greater 
contrast to sky in the near-infrared than in the optical. 

Infrared parallax observations have been pioneered by 
Jones (2000), who targeted the extremely active (and 
unfortunately at 76 pc also quite distant) late M dwarf 
PC0025+0447, as well as the nearby M dwarf VB10. The 
USNO also has an infrared astrometric program in oper- 
ation, from which published results are expected shortly 
(Vrba et al. 2002). 

3. OBSERVATIONS & SAMPLE 

Observations were carried out at 7 epochs over the pe- 
riod 2000 April 17 to 2002 May 30. At each epoch, observa- 
tions were carried out on either two half- or two full-nights. 
All observations were obtained with the SOFI infrared 
camera on the European Southern Observatory (ESO) 
3.5m New Technology Telescope (NTT). SOFI was used 
in its "large field" mode in which it provides a 4.92'x4.92' 
field-of-view with 0.28826" pixels (cf. Section 4.2). 

Exposures of each target were acquired with a fixed 
dither pattern (Fig. 2) as eight 120s exposures though 
the SOFI J filter. The exposure pattern was designed so 
that this 16 minutes of dithered exposure time sampled 
many different inter-pixel spacings. As much as was fea- 
sible (given observing time constraints) we attempted to 
acquire all epoch observations at the same hour angle as 
the very first epoch observation acquired, so as to minimise 
DCR effects. Each epoch observation was also carried out 
with a specified reference star positioned within a few pix- 
els of its position when observed on the very first epoch. 
This ensures all observations are carried out as near dif- 
ferentially as possible. 

Seeing conditions over the course of this program var- 
ied. Figure 1 shows a histogram of the seeing full- width at 
half maxmimum for all our astrometric observations. The 
median seeing was 0.82", with 80% of data being acquired 
in seeing conditions between 0.55 and 1.25". 

In addition to these epoch observations, all targets were 
also observed as they rose and set, so that DCR calibra- 
tions for each target could be developed, following the 
technique described in Tinney et al. (1995); Tinney (1996). 
The J filter was chosen for these observations as it of- 
fers the best contrast between sky- and T dwarf bright- 
ness. Typical near-infrared sky colours at La Silla are 
J-H=1.7, H-K=1.0-2.0 (2.0 for dark, 1.0 for bright) (ESO 
SOFI on-line documentation). By contrast typical T dwarf 
colours are -0.5<J-H<0.5 and H-K<0.5 (Burgasser et al. 



2002b), which means most T dwarfs are around a magni- 
tude brighter compared to the sky in J, than they are in 
H or K. 

The sample of objects observed is listed in Table 1, along 
with indications as to which targets were observed at which 
epochs. SOFI has a nominal gain of 5.6e-/ adu, and a nom- 
inal read-noise of 14e- per exposure. 

3.1. Object Names 

With the exception of elndB (Scholz et al. 
2003), all the objects discussed in this pa- 
per have been discovered by either the 2MASS 
(www.ipac.caltech.edu/2mass), SDSS (www.sdss.org) 
or DENIS (cdsweb.u-strasbg.fr/denis.html) sky sur- 
veys, and have been given object names by those surveys, 
based on their positions in J2000 coordinates. These 
names have the advantage of being very specific and infor- 
mative, and the disadvantage of being lengthy and clumsy. 
Throughout this paper, therefore, we will generally give an 
object's complete name when it is first used, and thereafter 
refer to it (when not confusing to do so) by a shortened 
2Mhhmm, SDhhmm or Dhhmm form where hh and mm 
are the right ascension hour and minute components of its 
name. 

4. ANALYSIS 

The analysis adopted for these data falls into two main 
areas: processing to produce linearised, flattened and sky 
subtracted images, which was quite specific to the SOFI 
instrument; and astrometric processing of these images, 
which identically follows that described in Tinney (1993); 
Tinney et al. (1995); Tinney (1996). 

4.1. Processing SOFI data 

Dark frames and zero-points : Dark frames obtained 
with SOFI reveal significant structure, which can be bro- 
ken down into a few components. 

1. a significant (^50-100 adu peak-to-peak) vertical 
structure, known as the "shade", which varies in 
intensity and shape with the overall level of illumi- 
nation of the array; 

2. a small (1-20 adu) dark current from the instrument 
and a small readout amplifier glow in each quadrant; 
and 

3. a tiny (<1 adu) but fixed "ray" pattern left after the 
previous two arc modeled and removed from dark 
current data. 

The shade pattern is of most concern as the remaining 
fixed patterns are small compared to the sky brightness. 
Figure 3 shows a dark current image displaying the shade 
effect, together with a set of vertical medianned profiles 
through the shade. Unfortunately, this shade profile is not 
constant - it varies with the intensity of the overall level 
of illumination of the array during an exposure, meaning 
one has an unknown zero-point for every pixel in every 
exposure. 

Calibration of the shade was achieved as follows. A 
small aperture (used to mask the instrument entrance 
when observing with one of the smaller fields of view) was 
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Table 1 

NTT Infrared Parallax Target Sample. 



Object 


Position 


Apr 


Jul 


Mar 


Apr 


Jul 


Mar 


May 


Ret 




(J2000) 


'00 


'00 


'01 


'01 


'01 


'02 


'02 




2MASS J05591914-1404488 


05 h 59 m 19? 1- 14° 04' 49" 


X 






X 




X 




a 


SDSS J10210969-0304201 


10 h 21 m 09?7-03°04'20" 




X 


X 


X 


X 


X 




b 


2MASS J10475385+2124234 


10 h 47 m 53?8+21°24'23" 


X 


X 


X 


X 




X 




c 


2MASS J12171110-0311131 


12 h 17 m ll?l-03°ll'13" 


X 


X 


X 


X 


X 


X 




c 


2MASS J12255432-2739466AB 


12 h 25 m 54?3-27°39'47" 


X 


X 


X 


X 


X 


X 




c d 


SDSS J12545390-0122474 


12 h 54 m 53?9-01°22'47" 




X 


X 




X 


X 




b 


SDSS J13464645-0031504 


13 h 46 m 46?4-00°31'50" 


X 


X 


X 


X 


X 


X 




e 


2MASS J15344984-2952274AB 


15 h 34 m 49 s 8 _ 29 o 52 '27" 


X 


X 


X 


X 


X 


X 




f d 


2MASS J15462718-3325111 


15 h 46 m 27!2-33°25'll" 


X 


X 


X 


X 


X 


X 




c 


SDSS J16241437+0029156 


16 h 24 m 14?4+00°29'16" 


X 


X 


X 


X 


X 


X 




g 



a Burgasser et al. (2000b) 
b Leggett et al. (2000) 
c Burgasser et al. (1999) 

d Shown to be a roughly equal binary by Burgasscr et al. (2003a) 
e Tsvetanov et al. (2000) 
•fBurgasser et al. (2002b) 
^Strauss et al. (1999) 
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Fig. 1. — Seeing histogram for all observations in this program. 
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Fig. 2. — Fixed dither pattern used for all epoch observations. Offsets arc in arcseconds from the initial base pointing position. 
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inserted into the NTT focal plane, and a series of flat fields 
obtained with varying exposure times. Because only the 
central quarter of the array is illuminated by this pro- 
cedure, it is possible to extract a shade profile from the 
edge of each image. It is also possible to record the level 
of illumination of the array which produced that shade 
profile. By performing a least-squares cubic polynomial 
fit through each pixel of these shade profiles (which cor- 
respond to rows on the detector) as a function of array 
illumination, it is possible to develop a parametrization 
for the shade profile. Using this parametrization it is a 
simple matter to produce a shade profile estimate for each 
data image, and subtract it. The result is an image with 
zero-point constant across the array. 6 

Linearity : All infrared detectors are non-linear to 
some extent. For SOFI, ESO usually recommends keeping 
sky and target object intensities below 10,000 adu in order 
to maintain linearity at better than 1%. Unfortunately, 
such an observing strategy is not useful for astrometry, 
which demands the largest possible dynamic range to en- 
sure targets (and reference stars) of widely differing mag- 
nitudes are usable in widely varying seeing conditions. We 
therefore calibrated the linearity of SOFI using the same 
shade profile data obtained above. Once the data have 
been shade corrected, they can then be used to exam- 
ine the response of each pixel to a constant light source 
over widely varying exposure times. Repetition of a 'cal- 
ibration' exposure time throughout the sequence allows 
the lamp's constancy to be calibrated - usually to within 
± 0.5%. A sample of the resulting linearity correction is 
shown in Fig. 5 for one of the SOFI quadrants. In all cases 
these tests were performed independently for each quad- 
rant, and the results were always consistent with the same 
linearity correction for all quadrants. A single correction 
was therefore derived as the mean of those in each quad- 
rant. Fig. 5 shows that the detector is «2.5% non-linear 
at 20,000adu above bias, but that data can be obtained 
and linearity corrected even up to 25,000adu. 

To linearise a pixel, then, with raw intensity Ty, it is 
simply necessary to multiply it by the polynomial P(JVj) 
= ciq + ai Iij + a>2 ifj + 0,3 ifj . The coefficients adopted were 
a = 1.0, ai = 0.0, a 2 = 1.11329 x 10- 10 ,a3 = -2.46799 x 
10~ 15 for 2000 April - 2001 April, and a = 1.0, ai = 
0.0, a 2 = 8.6124661 x 10- n ,a 3 = -1.6986849 x 10~ 15 for 
2001 July - 2002 May. 

Inter-quadrant Row Crosstalk : HgCdTe detec- 
tors typically show an effect known as inter-quadrant row 
crosstalk (Finger & Nicolini 1998). This has the effect that 
a constant, small fraction of the total flux seen in each row 
is seen as crosstalk at the same row in all the other quad- 
rants. Correction of this effect is straightforward. The de- 
tector is integrated up into a single vertical column, then 
the two halves of this cut (Y= 1-512 and Y=513-1024) are 
averaged, multiplied by a single cross-talk constant, and 
subtracted from every column of the detector. We found 
a crosstalk coefficient of 2.8xl0~ 5 worked well. 

Flat-fielding : Flat-fielding was performed using dome 
flats. Because of the (variable) shade pattern present in 



every dome flat, NTT staff have developed an observing 
recipe to obtain a "special" flat-field without a shade pat- 
tern present. Or one can used the shade calibration proce- 
dure described above to correct standard dome-flat fields. 
Both were tried for this program and both provided simi- 
lar results. In the end every run's data was flattened with 
a "special" dome flat, as per usual for SOFI observing. 

Sky subtraction : Each group of eight 120s dithered 
exposures was then used to create a normalised and me- 
dianned sky frame, which was re-normalised to each of 
the twelve observations to perform sky subtraction. So 
that the data frames would maintain approximate photon- 
counting errors, an appropriate constant sky level was then 
added back in to each frame 

Astrometric processing : Following this process- 
ing then, we have eight bias-subtracted, linearised, cross- 
talk corrected, flattened and sky-subtracted data frames 
for each astrometric epoch. These were then subject to 
further processing (ie. object finding and point-spread 
function fitting using DAOPHOT, DCR calibration and 
proper-motion and parallax solution fitting) as eight in- 
dividual observations, in a manner identical to that de- 
scribed in Tinney (1993); Tinney et al. (1995); Tinney 
(1996). 

4.2. Astrometric Calibration of SOFI 

Astrometric calibration observations were acquired in 
USNO Astrometric Calibration Region M (Stone et al. 
1999) on 2001 July 12 and 13. These consisted of six- 
teen 60s exposures (on each night) scattered throughout 
the 3.2° x7. 6° region which Stone et al. (1999) have as- 
trometrically calibrated. These were processed identically 
to our main astrometric targets. Reference catalogue po- 
sitions were extracted from the USNO ACR catalogue 7 
in SOFI-ficld-sized regions around each nominal telescope 
pointing position. These positions were then tangent pro- 
jected (using the SLALIB library Wallace 1999) to pro- 
vide reference data sets in arcsccond offsets on the sky 
for each observation. These were matched against the 
observed data to derive a set of linear (ie shift, scale 
and rotate) transformations from the SOFI pixel posi- 
tions to arcsecond offsets on the sky 8 . These transfor- 
mations determine the SOFI plate scale on this night to 
be 0.28826±0.00003"/pixel, and that the detector pixel's 
misalignment with N-S (0.030±0.003°). These data were 
also analysed to examine the amount of astrometric dis- 
tortion (ic. variability in the instrument plate scale with 
position in the field) present in the SOFI optical train. 
The astrometric calibration data show there is no signifi- 
cant astrometric distortion in the SOFI Large Field optics. 
The plate scale in the field corners is the same as that in 
the field center to within 0.1%. The SOFI field can be 
considered astrometrically flat to 0.1%. 

5. RESULTS FOR T DWARFS 

Astrometric solutions for our T dwarf targets were eval- 
uated in a manner identical to that used by Tinney et al. 
(1995). Briefly the procedure is to transform (using a lin- 



6 Sample paramctcrizations can be found at http://www.aao.gov.au/local/www/cgt/sofi. 

7 This data can be obtained from the Vizier service. 

8 This step made extensive use of M.Richmond's excellent match implementation of the Valdes et al. (1995) object list matching algorithm, 
which is available at http://acdl88a-005.rit.edu/match/ 
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(a) 





Fig. 3. — (a) Example of the SOFI shade pattern in a dark frame, (b) SOFI shade profiles obtained as vertical cuts through dark regions 
of the partially illuminated detector, for a range of count rates in the illuminated regions of the detector. The profiles span count rates in 
adu/pix from top to bottom of 3199, 2221, 1861, 1467, 1086, 958 701, 685, 554, 483, 226 and 87. 
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Fig. 4. — J-band finding charts with our target T dwarfs (cross hairs), and astrometric reference stars (boxes) highlighted. Each field is 
295.2" on a side. The direction bars in each field are 1' in size. 
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5000 10 1.5x10 4 2x10 4 2.5x1 4 

Observed odu above bias 

Fig. 5. — A linearity fit for a single SOFI quadrant. The data points plotted are the counts a perfectly linear detector would sec, divided 
by the observed counts, as a function of observed counts. The fit is a fourth order polynomial, which models the amount of non-linearity as 
a function of observed counts. The actual linearity correction coefficients used are given in the text. 





Fig. 6. — Plots of the observed astrometric motions of our target T dwarfs, relative to the ensemble of selected background reference stars. 
For each object we show plots of their motion in right ascension (Aoj) and dcclination(A<5) as a function of date. Also shown are the fitted 
proper motion + parallax (solid lines) and where a significant parallax has been measured, the fitted proper motion alone (dotted lines). The 
parameters fitted are given in Table 2. 
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Table 2 

NTT Infrared Parallax Sample - Relative Astrometry 



Object" 




7T r 


fir 


r 








(mas) 


(mas) 


( ) 


(km/s) 


2M0559 


71,5,11 


(96.9)" 


677.4±2.5 


122.6±0.1 


33.1±0.1 


S1021 


112,7,0 


O A A 1 A C 

34.4±4.6 


1 O O O 1 O A 

183.2±3.4 


njn O 1 1 A 

248.8±1.0 


25.2±2.4 


2M1047 


70,7,4 


110.8±6.6 


1698.9±2.5 


256.4±0.1 


72.7±4.4 


2M1217 


143,10,4 


90.8±2.2 


1057.1±1.7 


274.1±0.1 


55.2±1.4 


2M1225AB 


128,12,7 


75.1±2.5 


736.8±2.9 


148.5±0.1 


46.5±1.7 


S1254 


104,7,6 


73.2±1.9 


491.0±2.5 


284.7±0.1 


31.8±1.0 


2M1346 


118,9,5 


68.3±2.3 


516.0±3.3 


257.2±0.2 


35.8±1.4 


2M1534AB 


140,11,8 


73.6±1.2 


268.8±1.9 


159.1±0.1 


17.3±0.4 


2M1546 


150,10,9 


88.0±1.9 


225.4±2.2 


32.5±0.6 


12.1±0.4 


S1624 


152,11,8 


90.9±1.2 


373.0±1.6 


268.6±0.3 


19.5±0.3 



a See Table 1 for full object names 
6 Dahn ct al. (2002) 



ear transformation with rotation and a scale factor) all the 
frames for a given object, onto a chosen master frame of 
good seeing (known to have the detector rows and columns 
aligned with the cardinal directions to within ±0.1°), us- 
ing a set of well exposed reference stars which were re- 
quired to appear in every frame; differential colour refrac- 
tion (DCR - see Tinney et al. (1995)) coefficients were 
then evaluated for each of these reference stars (relative 
to the unknown 'mean' DCR coefficient for the reference 
star set), and the reference frame corrected for DCR; each 
frame was then re-transformed onto the 'master' frame; 
the DCR coefficient for the program object (relative to 
the DCR-corrected reference frame) was then evaluated; 
the program object was DCR corrected; and finally, an as- 
trometric solution in parallax and proper motion was made 
independently for both the a and S directions, using a lin- 
ear weighted-least-squares fit. Uncertainties arising from 
the DCR correction and the residuals about the reference 
frame transformation for each frame were carried through 
to this solution fit, so observations taken in poor seeing or 
with poor signal-to-noise due to cloud autmoatically re- 
ceive low weight. The final parallax was taken to be the 
weighted mean of the a and 6 solutions. Finding charts 
for our target stars (taken from our SOFI data) showing 
both the target and reference stars adopted can be seen in 
Fig. 4. 

The resulting relative astrometry is presented in Tabic 
2, the columns of which show; the number of frames (N/), 
nights (N„) and reference stars (N s ) used in each solution; 
the parallax and proper motion solutions (relative to the 
background reference stars chosen); and the derived tan- 
gential velocity for each target )based on the measured 
parallax, except for 2MASS0559, for which we adopt the 
distance of Dahn et al. 2002). Plots of these fits are shown 
in Fig. 6. The reference stars used to obtain this relative 
astrometry are typically within ±lmag. of the apparent 
magnitude of our target T dwarf. At these magnitudes 
(J=15-18) the reference stars will most commonly be G- 
to early M type stars at distances of 500-2000pc. Thus 
although we do not have the photometry available to es- 
timate detailed corrections from relative to absolute par- 
allax, we can estimate with some confidence, that such 



corrections will generally be less than 1 mas in size, and so 
not significant in comparison to our random astrometric 
uncertainties. 

It is instructive to examine the root-mean-square resid- 
uals obtained for the reference frame stars in our astrom- 
etry, since they tell us how precise we can expect the as- 
trometry of our target objects to be. Over the course of our 
program we found that for a single 120s exposure, the me- 
dian value of this rms residual was 0.042 pixel or 12.1 mas, 
with 80% of observations having an rms residual between 
6.9 and 20.2 mas. Recall that at each epoch we acquired 8 
such observations in a total exposure time of 960s, which 
would suggest the median precision from a single epoch is 
12.1/sqrt(8) = 4.3 mas. Residuals within these groups of 
eight were somewhat correlated (presumably because they 
are largely acquired in similar seeing conditions). 

The USNO have published astrometry for three T 
dwarfs: 2MASS0559, SDSS1254 & SDSS1624 (Dahn ct 
al. 2002). While all three were included in our pro- 
gram, insufficient epochs were obtained for 2MASS0559 
to measure a parallax. The equivalent relative parallax 
solution quantities for those we obtained (Table 2), are 
given by Dahn et al. as; for SDSS1254-01: 84.1±1.9mas, 
496.1±1.8mas/y, 285.2±0.4°, and for SDSS1624+00 : 
90.7±2.3mas, 383.2±1.9mas/y, 269.6±0.5°. These inde- 
pendent observations and solutions agree within uncertain- 
ties for almost all parameters - the exception being the 
parallax for SDSS1254, for which the two solutions are dif- 
ferent by about 5-ct, though Dahn et al. do comment that 
with only 1.2 years of data on this target their solution is 
only considered to be preliminary. 

Finally we note that with only 3 epochs of observation 
per year over two-plus years, there is always the possi- 
bility that systematic errors on individual runs may have 
impacted on our results. For example, a major change in 
SOFI's astrometric distortion or a decollimation between 
the telescope and SOFI on a single run could systemati- 
cally effect our results. The only way to detect such prob- 
lems is by detecting a poor match between our astrometric 
model and the data we obtain, which is difficult with less 
than 6 epochs. We believe the likelihood of this is small 
because; (1) exactly the same automated telescope image 
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analysis procedures were used to control the NTT's pri- 
mary figure throughout every night of every run, making 
the chance of an unusual NTT collimation with SOFI un- 
likely; (2) SOFFs astrometric distortion (as we have shown 
above) is tiny, so changes in it can have only negligible ef- 
fect; (3) SOFI is a Nasmyth mounted instrument, and so 
is always mounted horizontally and subject only to rota- 
tion about its optical axis, greatly reducing the likelihood 
of flexure within the instrument; and finally (4) because 
infrared instruments sit in temperature-controlled dewars 
they suffer almost none of the temperature-dependent flex- 
ure and defocus effects present in optical reimaging sys- 
tems, and they are also much less prone to being opened 
and modified over the course of an astrometric program. 
On-going monitoring and independent observations by in- 
dependent programs are the best way to test for unfore- 
seen systematic errors, and we look forward to checking 
our results against programs being carried out elsewhere. 

6. DISCUSSION 

6.1. Differential Colour Refraction for T dwarfs 

An interesting result of the DCR calibrations we per- 
formed for our T dwarf targets, was to find that T dwarfs 
have effective wavelengths in the J-band which are essen- 
tially indistinguishable from the ensemble of background 
reference stars against which their positions are mea- 
sured. This is shown in Fig. 7 which plots histograms 
of the DCR coefficients determined for reference stars and 
programme T dwarfs - the similarities in the ensemble 
values are clear. (These coefficients were derived using 
the method described in Tinney 1993 and Tinney 1996. 
Typical uncertainties in the individual determinations are 
w ±2 — 6mas/tan(ZA).) As a result, though we have cal- 
ibrated and applied DCR corrections to our data, such a 
procedure is not strictly necessary for near-infrared obser- 
vations of T dwarfs. These observations, therefore, are not 
rigidly tied to being carried out near the meridian, which 
adds enormously to the flexibility and efficiency of infrared 
parallax programmes. 

6.2. Photometry for L and T dwarfs 

There are currently only a few large and systematic pho- 
tometric databases for late M, L and T dwarfs extant. The 
first is the photometry from the 2MASS database Cutri et 
al. 2001, which has the advantage of being a well estab- 
lished photometric system which covers the whole sky, and 
includes all of our T dwarf targets, and almost of all of the 
other known L and T dwarfs (Burgasser 2002a). Unfor- 
tunately, the photometry for these objects in the J and 
K s 2MASS bands is often near the 2MASS photometric 
limits, so typical uncertainties of ±0.1 mag. or greater are 
not uncommon. Moreover, because 2MASS does not in- 
clude an optical passband, colour information has to rely 
on the J-K colour, which is typically small compared to 
the photometric precision, as well as giving only a small 
wavelength "lever arm" on the spectral evolution of L and 
T dwarfs. 

We make use of absolute Mj and Mk s values on the 
2MASS system compiled by Burgasser (2002a) which is 
based on the parallaxes presented in Dahn et al. (2002). 

9 Indeed the two arc so different that a distinctive distinctive name - 



Dahn et al. (2002) also present optical photometry in the 
Ic passband, and J,H,K photometry in a photometric sys- 
tem approximating that of the CIT system of Elias et al. 
(1982), as well as data form other work transformed onto 
this system. 

A second extensive database is that compiled by Leggett 
et al. (2002). This includes Z band photometry (on a 
UKIRT defined photometric system) as well as J,H,K pho- 
tometry transformed by the authors onto the MKO pho- 
tometric system (see Leggett et al. 2002, Section 3 for 
details). Because these data were acquired with a 4m tele- 
scope, their photometric precision is much higher than 
that for 2MASS. Great care should be taken in inter- 
comparing these two sets of photometry - the systematic 
differences between the two photometric systems are very 
significant. This is particularly true of the Z photometric 
system of Leggett et al. (2002), which is based on a rela- 
tively narrow interference filter (0.851-1.056 /im) used with 
a HgCdTc infrared array, leading to effective wavelengths 
for L and T dwarfs of «1.0/im, unlike the more common 
optical Z-type observations which are based on long-pass 
filters (>0.85^m) and the declining sensitivity of CCDs 
at >l^m, leading to effective wavelengths «0.9 /im. 9 . 
UKIRT Z photometry should not be assumed to be di- 
rectly comparable with optically based Z photometry In 
the discussion that follows, therefore, we will discuss only 
features in the absolute magnitudes within an individual 
photometric system. For this reason we do not make use 
of the more heterogeneous J,H,K compilation of Dahn et 
al. (2002). 

To these data sets we add observations of the re- 
cently announced T dwarf elndB Scholz et al. (2003), 
which has a mean I=16.7±0.1 and 2MASS photom- 
etry of J=11.91±0.04 and Ks=11.21±0.04 (Burgasser, 
priv.comm.). 

6.3. Photometric corrections for known binaries 

Several of the systems published in the photometric 
compilations listed above are known to be binaries, hav- 
ing been resolved either from the ground (Koerner et al. 
1999; Leggett et al. 2001), or using HST (Martm, Brander 
& Basri 1999; Reid et al. 2001; Burgasser et al. 2003a). 
Unfortunately, not all these systems have measured mag- 
nitude differences in all the passbands of interest, so we 
are forced to estimate magnitudes for the A and B com- 
ponents of these systems based on available colour-colour 
relationship data. In some cases (especially 2M0850B) 
these extrapolations are large, and the decomposed magni- 
tudes should be treated as indicative only. Table 3 shows 
the magnitude differences between each component and 
the total magnitude of each system, along with estimated 
spectral types from Burgasser et al. (2003a). Because 
of the similarity between the effective wavelengths of the 
HgCdTe-based Z and J bands, we assume that the mag- 
nitude differences in Z are the same as those derived at J. 
We do not differentiate here between UKIRT and 2MASS 
J and K bands. 

2M0746AB, 2M1146AB & 2M0850AB were observed by 
Reid et al. (2001) in the HST F814W filter from which 
magnitude differences for the components were derived in 
the Ic passbands. Infrared J-band magnitude differences 

- for these HgCdTe-based Z magnitudes is being widely adopted 
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□ T-dwarfs □ Reference Stars 



Fig. 7. — DCR coefficient histograms for reference stars and programme T dwarfs, in units of mas/Tan(Zenith Angle). See Section 6.1 for 
more details. 




Fig. 8. — Plot of Icr— J against J-Ks for L and T dwarfs from Burgasser (2002a) and Dahn et al. (2002). The data breaks into two clear 
sequences for the L and T dwarfs. The solid lines are least-squares linear fits to the data as follows: for L dwarfs J-Ks = 0.707(Ic-J) - 1.003 
for 3.2<I C -J<4.4 (<r = 0.14); for T dwarfs J-Ks = 7.598 - 1.371(I C -J) for 4.4<I C -J<6.0 (a = 0.25). 
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Table 3 

l and t dwarf binary magnitudes 



System 


Comp. 


Ale" 


AJ a 


AK a 


Sp.T.o 


Notes c 


2MASS J0746425+200032AB 


A 


0.50 


0.54 


0.59 


L0.5 


AIc,AJ R01. AK Fig. 8 




B 


1.12 


1.01 


0.95 


L0.5 


AIc,AJ R01. AK Fig. 8 


2MASS J1146345+223053AB 


A 


0.61 


0.64 


0.67 


L3 


AI C ,AJ R01. AK Fig.8 




B 


0.92 


0.87 


0.84 


L3 


AI C ,AJ R01. AK Fig.8 


2MASSs J0850359+105716AB 


A 


0.29 


0.55 


0.26 


L6 


AI C ,AJ R01. AK Fig.8 




B 


1.63 


0.99 


1.67 


T2: 


AI C R01. AJ Fig.10. AK Fig.8 


DENIS-P J0205290-115925AB 


A 


0.75 


0.75 


0.75 


L7 


Assumed equal mass binary K99,L01 




B 


0.75 


0.75 


0.75 


L7 


Assumed equal mass binary K99,L01 


DENIS-P J1228138-154711AB 


A 


0.54 


0.66 


0.75 


L5 


AJ M99, AK K99, AI C Fig.8 




B 


1.02 


0.86 


0.75 


L5 


AJ M99, AK K99, AI C Fig.8 


2MASS J12255432-2739466AB 


A 


0.24 


0.28 


0.16 


T6 


AI C ,AJ B03, AK Fig.8 




B 


1.76 


1.63 


2.15 


T8 


AI C ,AJ B03, AK Fig.8 


2MASS J15344984-2952274AB 


A 


0.53 


0.75 


0.75 


T5.5 


AI C ,AJ B03, AK assumed 0.75 




B 


1.03 


0.75 


0.75 


T5.5 


AI C B03, AJ,AK assumed 0.75 



a Difference in magnitude between the component and the total magnitude of the system in this passband. 
b Burgasscr et al. 2003a 

c R01 - Reid et al. 2001, K99 - Koerner ct al. 1999, L01 - Lcggctt ct al. 2001, M99 - Martin, Brander & Basri 1999, B03 - Burgasser et al. 



were then estimated using the L dwarf sequence of the 
Mj versus Ic~J colour-magnitude diagram (which has a 
roughly constant slope). With the exception of 2M0850B, 
this procedure will be adequate for all the L dwarfs, and 
those values are shown in Table 3. 2M0850B is an excep- 
tion because its absolute magnitude at Iq is so faint that it 
must be an early- to mid-T dwarf, rather than an L dwarf. 
And, as we show in Section 6.5, the colour-magnitude di- 
agram is not even remotely linear across the L-T dwarf 
transition. For 2M0850B, therefore, we have used the Mj 
for the AB system of Dahn ct al. (2002), and the magni- 
tude differences of Reid et al. (2001) to derive for the B 
component M/ = 20.02±0.23. The colour-magnitude di- 
agrams in Section 6.5 then imply Ic~J ~ 5.1±0.2, from 
which we derive the B component J magnitude difference 
shown in the Table. 

To derive K-band magnitudes for the components of 
these systems, we have plotted J-Ks versus Ic~J for all 
the L and T dwarfs in Burgasser (2002a) and Dahn et al. 
(2002) in Figure 8. The data reveal two separate sequences 
- the L dwarfs in which J-Ks becomes redder along with 
Ic-J, and the T dwarfs in which the reverse holds. The 
two lines on the plot are linear fits to these two regimes 
(arbitrarily divided at ^^=4.4). From these relations 
we predict J -Ks colours for L and T dwarfs from their I- J 
colours, and so derive the magnitude differences for each 
component in Table 3. 

D0205AB & D1228AB were both observed by Koerner 
et al. (1999) in the K band at Keck, and D0205AB was in- 
dependently observed at UKIRT by Leggett et al. (2001). 
D1228 was also observed in the J band with HST by 
Martm, Brander & Basri (1999). D0205 was found to be 
a pair of objects with equal brightness at K, and in the 
absence of any other information we assume it to be an 
equal mass binary. D1228 is a nearly equal mass binary - 
from the marginal J-K colour difference between the two 
components we can extrapolate to a magnitude difference 
between the components at I of 0.48. 

2M1225AB, 2M1534AB have been observed by Bur- 
gasser et al. (2003a) with HST in the F814W and F1042M 
filters (the latter enabling the derivation of approximate 



J magnitude differences for the systems). Once again we 
use the I-J colours of these objects to extrapolate to K 
magnitude differences for 2M1225AB's components. For 
2M1534AB the magnitude differences estimated at F1042 
are only marginally different from zero, so we assume equal 
brightnesses in this system at J and K. 

6.4. Spectral-Type-Magnitude relations for L and T 

dwarfs 

Table 4 lists Burgasser (2002a), Dahn et al. (2002) and 
Leggett et al. (2002) photometry for our NTT parallax 
sample, along with the resulting absolute magnitudes in 
these systems. Also listed are spectral types on the scheme 
of Burgasser et al. (2002b). 

Figure 9 shows plots of spectral type against Mz, Mj, 
Mj and M^ s /M^. Also shown are absolute magnitudes 
for late M and L dwarfs using parallaxes and 2MASS pho- 
tometry from Dahn et al. (2002) for the 2MASS panels, 
and parallaxes from Dahn et al. (2002) and UKIRT pho- 
tometry from Leggett et al. (2002) for the UKIRT panels. 
The spectral types are on the system of Kirkpatrick et 
al. (1999) for the M and L dwarfs, and Burgasser et al. 
(2002b) for the T dwarfs. Known multiple systems are 
noted with circles, and decomposed into their component 
magnitudes as discussed above. 

The two K-band plots (Fig. 9a and 9b) indicate that 
in both systems, the L-T transition is marked by a steep- 
ening of the spectral-type-magnitude relation. In general, 
however, the relationship between absolute magnitude at 
K and spectral type is well behaved for the purpose of es- 
timating absolute magnitudes from spectral types. This 
is certainly not true in the 2MASS and UKIRT J pass- 
bands (Fig. 9c and 9d). Indeed, both sets of data in- 
dicate a strong inflection (a "hump") in the relationship 
between absolute magnitude and spectral type for early 
T dwarfs - as a class, the T0-T4 brown dwarfs have ab- 
solute magnitudes brighter than the latest L dwarfs by a 
magnitude or more. Put another way, a simple extrap- 
olation of the spectral-type-magnitude relationship for L 
dwarfs (eg. that from Dahn et al. (2002) shown in the fig- 
ure) underestimates the absolute magnitude of the early- 
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Fig. 9. — Absolute magnitudes in a variety of bandpasses, as a function of L and T spectral types on the Kirkpatrick ct al. and Burgasser 
et al. systems. Panels (a) and (c) show absolute Mj and M^ s from 2MASS photometry in Burgasser (2002a); Scholz ct al. (2003) and 
this paper. Panels (b),(d) and (f) show absolute Mj, Mj( and Mz on the Leggett et al. (2002) UKIRT system. Panel (c) shows M; c from 
photometry in Dahn et al. (2002). All parallaxes are from either this work (open squares) or Dahn ct al. (2002) & Scholz ct al. (2003) (filled 
squares). The dotted line in panels (c) and (d) is the spectral- type-magnitude relation of Dahn et al. (2002). Known unresolved binaries are 
encircled, and decomposed into their resolved components as described in Section 6.3. The solid heavy lines are seventh-order polynomial 
weighted fits to the data (see Section 6.8 for coefficients). The decomposed JK magnitudes and spectral type of 2M0850B arc much more 
uncertain than the formal uncertainties plotted would indicate, and should be treated as indicative only. 
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Table 4 

T dwarf Photometry & Absolute Magnitudes for NTT Parallax Sample 



Object 


T" 


2MASS" 




MKO c 




2MASS 




MKO 








J 


K 3 


J 


K 


Z-J 


Mj 


M Ks 


M z 


M,, 


M K 


SO1021 


T3 


16.26±0.10 


15.10±0.18 


15.88 


15.26 


1.78 


13.94±0.29 


12.78±0.33 


15.34±0.27 


13.56±0.27 


12.94±0.27 


2M1047 


T6.5 


15.82±0.06 


16.30±0.30: 


15.46 


16.20 


1.93 


16.05±0.14 


16.52±0.33 


17.61±0.13 


15.68±0.13 


16.42±0.13 


2M1217 


T7.5 


15.85±0.07 


15.90±0.30: 


15.56 


15.92 


2.00 


15.64±0.09 


15.69±0.30 


17.35±0.06 


15.35±0.06 


15.71±0.06 


2M1225 


T6 


15.22±0.05 


15.06±0.15 


14.88 


15.28 


1.89 


14.60±0.09 


14.44±0.17 


16.15±0.08 


14.26±0.08 


14.66±0.08 


SD1254 


T2 


14.88±0.04 


13.83±0.06 


14.66 


13.84 


1.74 


14.20±0.07 


13.15±0.08 


15.72±0.06 


13.98±0.06 


13.16±0.06 


2M1346 


T6 


15.86±0.08 


15.80±0.30: 


15.49 


15.73 


2.24 


15.03±0.11 


14.97±0.31 


16.90±0.08 


14.66±0.09 


14.90±0.08 


2M1534 


T5.5 


14.90±0.04 


14.86±0.11 


14.60 


14.91 




14.24±0.05 


14.19±0.12 




13.94±0.05 


14.25±0.05 


2M1546 


T5.5 


15.60±0.05 


15.42±0.17 








15.32±0.07 


15.14±0.18 








SD1624 


T6 


15.49±0.06 


15.40±0.30: 


15.20 


15.61 


2.12 


15.28±0.07 


15.19±0.30 


17.11±0.04 


14.99±0.04 


15.40±0.06 



a Spectral types on the Burgasscr ct al. (2000a) system 

6 Burgasscr (2002a). K s entries marked ":" are upper limits. 

c Leggett et al. (2002), except for 2M1534 which is Leggett, priv.comm. J,K uncertainties typically ±0.03, Z-J typically ±0.05 



to mid-T dwarfs by up to two magnitudes. This "early 
T hump" has been noted previously (Dahn et al. 2002), 
though on the basis on fewer T dwarf parallaxes. It has 
been suggested (Burgasser 2002a) that binarity could be 
the cause for early T dwarfs being more luminous than the 
late L dwarfs. While it is certainly true that the L and T 
dwarfs which have been resolved as binaries are displaced 
to apparently high absolute magnitudes when plotted as 
unresolved objects, the addition of new parallaxes would 
seem to indicate the over-luminosity of early T dwarfs is 
a general property, rather than being due to the selection 
of objects which happen to be binaries. Moreover, the 
magnitude or more of over-luminosity is too large an ef- 
fect to be due equal-mass binarity which can produce a 
brightening of only 0.75 mag. A similar (though possibly 
less pronounced) inflection is seen in the Mz relation (Fig. 
9f), while the Mj relation (Fig. 9e) would appear to be 
almost as monotonic as that at K, though with a more pro- 
nounced inflection at the L-T boundary. Having said this, 
however, 2M0559 continues to appear to be over-luminous 
compared to the other early- to mid-T dwarfs in the figure. 
Burgasser et al. (2003a) failed to resolve a binary compan- 
ion in this system with HST, implying that if it is a binary 
it must have a separation of less than 0.5 a.u. We also note 
that it has been suggested (Tsuji & Nakajima 2003) that 
the selection of preferentially young objects could produce 
the "early T hump" - we discuss this further in Section 
6.6.2. 

There are good physical reasons for expecting a mono- 
tonic relationship between effective temperature (T c //) 
and luminosity (L) in these objects, since these quanti- 
ties are directly determined by interior (rather than pho- 
tosphcric) properties. However, it must be remembered 
that as proxies for T c // and luminosity, absolute magni- 
tude in a given passband and spectral type are far from 
perfect. Spectral typing is in essence an arbitrary alloca- 
tion of a quantity to an object based on what its spectra 
look like - there is no guarantee that the relationship be- 
tween spectral type and T e ff (even if monotonic) should 
not have significant changes in slope. Similarly, the rela- 
tionship between absolute magnitude in a given passband 
and luminosity is even more problematic. From the spec- 
tra of objects ranging from L to T spectral types, and in- 
deed from their J-K colours (Dahn et al. 2002), we know 
that significant changes take place in their photospheres. 
There is significant redistribution of flux in the spectra of 



brown dwarfs across the L-T transition. We should not 
be surprised if this results in the relationship between lu- 
minosity absolute magnitude in a given passband not only 
containing changes in slope, but not even being monotonic. 

Given our current parallax database, spectral type is a 
very poor proxy for absolute magnitude in the Z and J 
bands from mid-L to mid-T spectral types. The sequences 
in Fig. 9 will need to be filled in by many more L and 
T dwarfs before precise absolute magnitudes can be esti- 
mated from spectral types with confidence. 

6.5. Colour Magnitude Diagrams for L and T dwarfs 

Using the same photometry, we can construct a va- 
riety of colour-magnitude diagrams. Figure 10 shows 
such diagrams based around Cousins I c , UKIRT Z and 
both UKIRT and 2MASS J,K photometry, while Fig- 
ure 11 shows similar diagrams for UKIRT and 2MASS 
J-K colours. The most noticeable feature of these dia- 
grams is how few are actually useful as traditional colour- 
magnitude diagrams - almost none show the simple mono- 
tonic relationships between absolute magnitude and colour 
which hold for stars and brown dwarfs down to the early 
L dwarfs. 

Fig 10c, d show that I-K colours jumps to the blue by 
I-Kw0.5 mag as the L-T transition is crossed at M/«19, 
Mx~ 13, but then tends redward again for later and later 
T dwarfs. However, G1570D, one of the latest and faintest 
T dwarfs currently known, never becomes as red as the 
latest L dwarfs. This blueward jump is particularly pro- 
nounced at M/ where the absolute magnitudes of L8 and 
early T dwarfs arc indistinguishable. As a result I-K 
should be considered a poor indicator for determining the 
absolute magnitude or effective temperature of late L to 
late T dwarfs. In particular, and luminosity function based 
on I-K>5 will be subject to serious biases which will in- 
troduce completely spurious structure into the luminosity 
function. 

Fig 10a, b shows that I-J colour-magnitude diagrams can 
be considered the "best of a bad bunch" when it comes to 
the traditional use of colour-magnitude diagrams (i.e. es- 
timating absolute magnitudes from photometric colours), 
since the cooling curves of brown dwarfs do not reverse in 
I-J as they do for every other panel of Figs 10 and 11. 
Even so, between I-J=4 and I-J=5 they show the same 
pronounced "S-curve" seen in the spectral type data of Fig. 
9, with early T dwarfs being up to a magnitude brighter 
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Fig. 10. — Colour magnitude diagrams based on lc and 2MASS J,Ks photometry (Dahn et al. 2002; Burgasser 2002a)), and the parallaxes 
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in Mj than late L dwarfs. And, once again we see that 
2M0559 appears anomalously bright, suggestion binarity 
in spite of Burgasser et al. (2003a) 's failure to resolve it 
with HST.. 

The Z-J colour-magnitude diagrams (Fig. 10c, f) reveal a 
very steep colour-magnitude relation, with scatter which 
is significantly larger than the photometric errors. The 
slope of the colour-magnitude relation is so steep that no 
meaningful estimate of Mz or Mj can be derived from a 
Z-J colour. This is not surprising, given the very close 
effective wavelengths of HgCdTe-based Z and J photom- 
etry. There is some evidence for trend at the bottom of 
this colour-magnitude diagram that at Mz > 16.5 and 
Mifs ^ 14.5, Z-J colours becomes bluer for fainter and 
later- type objects. 

Colour magnitude diagrams involving Z-K (Figs lOg-h) 
and J-K (Figs lla-d) show an especially pronounced rever- 
sal of the brown dwarf cooling curves beyond Mr- « 12.5, 
Mj w 14 and M z * 16. This has been noted in J-K 
colour-magnitude diagrams by several authors (e.g. Bur- 
gasser et al. 2002c). Major changes take place in photo- 
spheres below the L-T transition, with the result that T 
dwarfs swap from very red, to very blue, J-K colours. It is 
interesting that the Z-K diagrams show almost identical 
behavior, with Z-K colours for the very faintest T dwarfs 
becoming as blue as Z-K«-l. This compares with colours 
based on the SDSS z' filter (see eg.Dahn et al. 2002 Fig. 3) 
which continue to become redder for the latest T dwarfs. 
This once again clearly demonstrates the considerable dif- 
ference between CCD-based z', and the HgCdTe-based Z. 

There is a clear warning to astronomers implicit in 
these diagrams - conclusions reached about luminosity- 
and mass- functions based on luminosity and/or colours for 
the L-T effective temperature range are fraught with diffi- 
culty. In particular, luminosity functions determined from 
the colours of objects in field samples will produce com- 
pletely spurious features in the derived luminosity- and 
mass-functions, unless the various "bumps and wiggles" in 
these diagrams are adequately and correctly modeled. (See 
for example Reid & Gizis (1997)'s demonstration of the 
formation of a "false" peak in M dwarf luminosity function 
based on the traditional - and inadequate - parametriza- 
tion of the M dwarf colour-magnitude relation) . Similarly, 
determining bolometric luminosity functions from appar- 
ent magnitudes in cluster-based samples is problematic, as 
we can expect similar "bumps and wiggles" to be present in 
the bolometric correction relations for the L and T dwarfs. 

Features like these can introduce significant systematic 
biases into the mass functions derived from even a per- 
fect statistical sample. Actual statistical data with all the 
added complexities of uncertain age and binarity distribu- 
tions add yet more complications. Monte Carlo simula- 
tions are essential to the interpretation of any luminosity- 
or mass-function in the L-T effective temperature range, 
it is important to carefully "reverse" model such functions 
from sets of mass-function models, through a variety of 
possible colour-magnitude and bolometric-correction rela- 
tions (as allowed by the extant data), to sample observa- 
tional data. Then such artificial data can then be meaning- 
fully compared to statistical samples in the observational 
plane. Mass- or luminosity functions which do not include 
such extensive reverse modeling should be treated with the 



utmost suspicion. 

6.6. Theoretical Models for L and T dwarfs 

Ultra-cool dwarfs are notoriously difficult to model - 
the components which need to be included in models for 
L and T dwarfs include (Allard et al. 1997): the effects of 
tens of millions of molecular transitions in species includ- 
ing H 2 0, CH 4 , TiO, VO, CrH, FeH, and a host of others; 
complex treatments of the line wings of enormously H 2 
and He pressure-broadened neutral alkali lines like K and 
Na, Rb and Cs; collision-induced molecular H2 opacity; 
both the chemistry and opacity involved in the conden- 
sation, settling, revapourisation and diffusion of a variety 
of condensates including liquid Fe, solid VO, and a range 
of aluminium, calcium, magnesium and titanium bearing 
refractories; and finally (and least readily modeled of all) 
the effects of rotation-induced weather on the cloud decks 
which condensates will form. 

Significant progress has been made in recent years on the 
detailed solution of photospheric models using very large 
line lists (see Allard et al. (1997) for a review). Probably 
the largest outstanding problem for modelers of L and T 
dwarfs is dealing with condensation. Three approxima- 
tions to this complex situation have currently been im- 
plemented. "Dusty" models (eg. the DUSTY model of 
Allard et al. 2001) assume condensates remain well sus- 
pended and in chemical equilibrium where they form in 
the photosphere. In general such models have been shown 
to work reasonably well for L dwarfs, suggesting that their 
cloud layers lie within their photospheres. "Condensation" 
models (eg. the COND models of Baraffc et al. 2003, and 
the CLEAR models of Burrows et al. 1997) neglect dust 
opacities, to simulate the removal of all condensates from 
the photosphere as they form (most likely through gravi- 
tational settling). The "CLOUDY" models of Ackerman 
& Marley (2001) and Marley et al. (2002) incorporate a 
model for condensate cloud formation, based on an as- 
sumed sedimentation efficiency parameter f ra m- 

6.6.1. Colour-magnitude diagrams in J-K and Z-K 

Figures 10 and 11 have over-plotted on them a variety of 
these models, including the DUSTY (Chabrier et al. 2000; 
Allard et al. 2001) and COND models (Baraffe et al. 2003) 
for an age of 1 Gyr, and the CLEAR and CLOUDY mod- 
els as presented in Burgasser et al. (2002c) for f ra in = 3 
(determined as the best fit for this model in Jupiter's am- 
monia cloud deck Marley et al. 2002). As previous studies 
have shown, DUSTY models reproduce the general fea- 
tures (if not the precise colours) of the cooling curves for 
L dwarfs, but then proceed to much redder colours than 
are observed beyond L8. This has been interpreted as in- 
dicating that condensates are present in the photosphere 
of L dwarfs. The COND and CLEAR models reproduce 
the general features of the cooling curves for late T dwarfs, 
indicating that at these effective temperatures condensate 
opacities do not contribute to the radiative transfer, which 
suggests that the condensate layers have dropped below 
the photosphere. DUSTY and CLEAR/COND models, 
therefore, describe the "boundary conditions" to the con- 
densate opacity problem, and are appropriate for the L and 
late T types respectively. But what about the intermediate 
case which must be appropriate to early T dwarfs? This 
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is exactly the situation with the sedimentation models of 
Marley et al. (2002) should be able to address. 

Burgasser ct al. (2002c) compared their CLOUDY 
models for f ra in — 3 with a 2MASS Mj:J-K colour- 
magnitude diagram (as we do in Fig. 11). As for 
the DUSTY models, the CLOUDY models predict the 
general behaviour of L dwarfs, and then veer towards 
bluer J-K colours at late T dwarf temperatures. How- 
ever, this transition does not match the observed se- 
quence, which transitions nearly horizontally between 
the L dwarf/DUSTY/CLOUDY sequence and late-T 
dwarf/CLEAR/COND model at Mj«14, M Ks ^13. (We 
note that though the equivalent models are not available 
in the UKIRT Z,J,K bandpasses, very similar behaviour is 
seen in Fig. 10, with a clear transition between the L dwarf 
and late T dwarf sequences.) Burgasser et al. (2002c) 
suggested that a possible resolution for this discrepancy 
could be the appearance of uneven cloud cover on the sur- 
face of early T dwarfs. This would allow the emergent 
spectrum to appear as a "mixture" of the CLEAR/COND 
and CLOUDY spectra. They modeled this by interpolat- 
ing between their CLEAR and CLOUDY models at effec- 
tive temperatures of 800 K, 1000 K, 1200 K, 1400 K, 1600 
and 1800 K with varying fractions of the two models (ic. 
20%,40%,60 & 80%). The tracks for these "mixture" mod- 
els are shown in Fig. 11 as dotted lines, and suggest that 
there is a transition sequence between L and T dwarfs at 
T c// « 1300K. SDSS1021, SDSS1254, 2M1225, elndB and 
possibly 2M0850B (though with some uncertainty because 
of the poor quality of its decomposed secondary flux) fill 
out this transition region. The status of 2M0559 is un- 
clear. If it is a single object, then it probably represents 
the 'top' of the late T dwarf cooling sequence, which is 
>1 mag. brighter in Mj than the bottom of the L dwarf 
sequence. If however, it is a binary, then the prototype for 
the 'top' of the late T dwarf cooling sequence is probably 
more like the object 2M1225A or 2M1346 at a spectral 
type of T5.5-T6. The "transition temperature" indicated 
by the additional T dwarfs in this work is slightly warmer 
(«1300K) than that found by Burgasser et al. (2002c). 

An alternative dust model to the sedimentation models 
of Marley et al. has been developed by Tsuji & Naka- 
jima (2003). These "Unified Cloudy Models" (UCM) are 
built around a single thin dust layer in which particles 
of size greater than a critical radius are removed from 
the photosphere by sedimentation. This critical radius is 
parametrized by a critical temperature T cr below which 
dust particles sediment, which is determined by comparing 
model results to colour-magnitude diagrams. This single 
model has the advantage of predicting the gross behaviour 
of brown dwarfs as they transition from L to T spectral 
types, with a single model. Unfortunately, (Tsuji & Naka- 
jima 2003, Fig. 2) the detailed behaviour of the models 
does not match observations. In particular, UCM cannot 
make L dwarf as red or faint in Mj:J-K as they actually 
appear. Nor does it predict the observed brightening of 
the "early T dwarf hump" other than as an age-selection 
effect, which we conclude below, is not the case. 

It should be noted, however, that the interpretation of 
Marley et al. models and data in Fig. 11 in terms of 
cloud openings (i.e. as providing evidence for the exis- 
tence of weather in early T dwarfs) is quite dependent 



on the details present in the Marley et al. (2002) models. 
An independent test of this conclusion is clearly desirable. 
Fortunately, Fig. 11 indicates that J and K band time- 
series photometry can provide that test. The location of 
a given object on the "transition sequence" will depend 
critically on its fractional cloud cover. Because this could 
be expected to change as each brown dwarf rotates, a sta- 
tistical study of the J-band variability from late L dwarfs 
to late T dwarfs should find stronger variability in early T 
dwarfs, than in late L dwarfs or late T dwarfs. 

Finally, we note that although the COND models do 
not do a very good job of predicting the absolute colours 
of late T dwarfs in Z-J and Z-K (Fig. lOc-h), they do 
suggest a trend for late T dwarfs to become bluer in both 
Z-J and Z-K as they get colder and fainter than Mj«14.5 
and M#-«14.75. Moreover, the available data suggest this 
trend is real, though the absolute colours of T dwarfs at 
these magnitudes are somewhat redder than the models 
would predict. 

6.6.2. The "Early T Hump" 

The Mj:I<7-J colour-magnitude diagrams shown in Fig. 
lOa-b indicate a remarkable brightening at Mj for the ob- 
served early T dwarfs. Unfortunately, neither the COND 
nor the DUSTY models indicate why this should be so. 
The DUSTY models predict an extension of the L dwarf 
sequence, which we have good reason to believe is not cor- 
rect, based on the analysis of colour-magnitude diagrams 
in J-K above. Unfortunately, the COND models also fail 
to look even remotely like the available data for T dwarfs in 
Fig. lOa-d. Shortcomings in these models at short wave- 
lengths have been noted by Baraffe et al. (2003), which 
are thought to be due to an inadequate treatment of the 
extremely broad wings of the K and Na lines at these wave- 
lengths. 

One possible interpretation of the "early T hump" in 
Fig. lOa-b, is that it could be a gravity effect (Tsuji 
& Nakajima 2003). Very young brown dwarfs will have 
isochrones slightly offset to brighter magnitudes than older 
brown dwarfs, because of their lower gravities. This effect 
is particularly pronounced in photospheres in which dust 
is an important opacity source. It is possible then that 
the "early T hump" could be produced by the preferential 
selection of young, bright brown dwarfs. 

Fig. 12 plots the same data as that shown in Fig. 10a- 
b, but now we plot four isochrones spanning 50Myr-5 Gyr 
to examine the effects of age. The figure shows that, 
as expected, the DUSTY models (most appropriate for 
L dwarfs) show significant offsets in their isochrones of 
a magnitude or more between 50Myr and 5 Gyr. These 
offsets are not as marked for the COND models. Unfortu- 
nately, interpreting the "early T hump" as an age effect is 
severely complicated by the fact that it occurs at exactly 
the point where there is good evidence to believe neither 
the DUSTY or COND models are working. 

For the L dwarfs and late T dwarfs, the spread in the 
colour-magnitude diagram is not pronounced, (particu- 
larly when known binaries are decomposed), suggestive of 
the small 100 Myr - 1 Gyr age spread seen in other studies 
of L and T dwarfs (Dahn et al. 2002; Scholz et al. 2003). 
It is certainly nowhere near as pronounced as would be 
necessary to produce the age spread required to account 
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Fig. 12. — Colour magnitude diagrams based on Iq and 2MASS J photometry (Dahn et al. 2002; Burgasscr 2002a)), and the parallaxes of 
this paper and Dahn et al. (2002); Scholz et al. (2003). DUSTY (Chabrier et al. 2000; Allard et al. 2001) (rightmost lines) and COND Baraffe 
et al. (2003) (leftmost lines) are shown as in Fig. 10, but now for a range of ages: 50 Myr (solid), 100 Myr (dot- dashed), 1 Gyr (dotted), 5 Gyr 
(dashed). 



for the more than 1 magnitude brightening of the "early 
T hump" all by itself. Moreover, there is a definite spec- 
tral type trend along the track represented by the "early T 
hump" , as seen in the spectral-type-magnitude diagrams of 
Fig. 9 - from the late L dwarfs, through elndB, SD1254, 
SD1021 to 2M0559. This same trend is seen in the colour- 
magnitude diagrams. We interpret this as indicating that 
the "early T hump" truly is a feature in the cooling curve 
of brown dwarfs, rather than an artifact of youth and se- 
lection. 

6.7. The Onset of CH4 absorption in Clusters 

Methane filters centered on the strong CH4 absorption 
bands in the H-band have been acquired by a number of 
observatories for use in their infrared cameras. Given we 
have now measured just where, in absolute magnitude, the 
T dwarf class occurs, the question arises, "At what mag- 
nitudes will CH4 absorption in young star clusters set in?" 
Fig. 13 shows UKIRT Mj and M H :3-K colour-magnitude 
diagrams, along with the DUSTY and COND models at 
ages from 10 Myr to 5 Gyr. Based on these diagrams, we 
can conclude that for field T dwarfs, as discovered by the 
2MASS and SDSS surveys, CH4 absorption (correspond- 
ing to spectral classes around T2 and later) sets in at 
at Mjy«13, and somewhat more confusingly at Mj~14 

though because of the brightening of "early T hump" 
at J non-CH4-absorbing L dwarfs will actually be fainter 
than the earliest CH4 absorbing T dwarfs. 

Because the turn on of CH4 absorption is primarily an 
effect driven by effective temperature, to first order it will 



occur at the same absolute magnitude in young clusters as 
it does in the field. Looked at in slightly more detail, how- 
ever, we can see that for a given colour in Fig. 13, there is 
a small offset to brighter magnitudes for younger objects 
- in the DUSTY atmosphere case a 10 Myr dwarf at the 
end of the L8 sequence will be fsl.Omag. brighter than a 
1 Gyr dwarf of the same colour and effective temperature, 
and 1.3 mag brighter than a 5 Gyr dwarf. In the COND 
case the equivalent differences are 1.3 and 1.7 magnitudes. 
The likely ages for our field T dwarfs will be somewhere in 
the range lOOMyr-lGyr (Dahn et al. 2002; Scholz et al. 
2003). This would suggest that in clusters like IC 2391 or 
IC2602 of age 10-20 Myr at d«150pc the absolute mag- 
nitude for CH4 onset will be M#~12-12.5, or equivalcntly 
H~18-18.5. For older clusters like the Pleiades (100 Myr, 
d^T25pc) these numbers are more like M#~12.5-13 or 
H^18-18.5.Both of these are eminently reachable magni- 
tude limits with wide-field cameras on 4m-class telescopes, 
suggesting that CH4 imaging may be a powerful tool for 
easily conducting an unbiased census of T dwarfs in large 
open clusters. Similarly for more compact, but distant 
clusters, like Trapezium (~25Myr, 450 pc) observations 
at 14^20-20.5 are tractable over the fields-of-view required 
on 8m-class telescopes. 

Fig. 13 also has implications for the interpretation of 
potential cluster membership. For example, Zapatero Os- 
orio et al. (2002) have found a T6 dwarf in the direction of 
the o Orionis cluster. Fig. 13 suggests that a field brown 
dwarf of this spectral type will have M#=15.0±0.5. For 
the much younger age of a Orionis (1-8 Myr Zapatero Oso- 
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Fig. 13.— Colour magnitude diagrams in UKIRT Mj and M H vs J-K. DUSTY (Chabrier et al. 2000; Allard et al. 2001) (rightmost lines) 
and COND Baraffc ct al. (2003) (leftmost lines) are shown as in Fig. 10, but now for a range of ages: 10 Myr, (heavy solid) 50Myr (solid), 
100 Myr (dot- dashed), 1 Gyr (dotted), 5 Gyr (dashed). The decomposed J and H magnitudes of 2M0850B are much more uncertain than the 
formal uncertainties plotted would indicate, and should be treated as indicative only. 



rio et al. 2002) this will be more like M ff =14.0±0.5, which 
would imply a distance to the a Ori J053810. 1-023626 T6 
dwarf of 192±50pc - more consistent with being a fore- 
ground object than a member of the cluster at d=352 pc 
(Perryman et al. 1997). 

6.8. Colour-Magnitude & Spectral-Type-Magnitude 
relations for L and T dwarfs 

Figures 9 and 10 have over-plotted on some of their pan- 
els high-order polynomial fits to the weighted (and binary 
decomposed) data. As inspection of the figures shows, 
these fits are not always particularly successful at mod- 
eling the extremely complex behaviour of these cooling 
curves in the observed passbands. Nonetheless, in the ab- 
sence of working atmospheric models the fits may be a 
useful tool, so long as their weaknesses are acknowledged. 
We therefore provide the coefficients for these fits, and the 
root-mean-square scatter about the fits in Table 5. 

7. CONCLUSION 

We have shown that high precision parallaxes can be 
obtained with common-user near-infrared cameras using 
techniques very similar to those used in optical CCD as- 
trometry. The new generation of large format infrared 
imagers based on HAWAII1 (IK) and HAWAII2 (2K) 
HgCdTe arrays, and the new generation of large format 
InSb arrays, offer exciting prospects for the astromctry of 
cool brown dwarfs in the future. Due to their infrared 
methane absorption bands, T dwarfs have quite similar 
effective wavelengths to the ensemble of background refer- 
ence stars, which makes the correction of differential colour 



refraction effects considerably easier. The "early T hump" 
(i.e. the brightening at the J band of early T dwarfs rel- 
ative to late L dwarfs) appears to be a feature of brown 
dwarf cooling curves, rather than an effect of binarity or 
age. And finally, these data imply that detection of T 
dwarfs in clusters can be made directly at tractable mag- 
nitudes in the H-band, opening the way to a new genera- 
tion of cluster mass function studies based on the powerful 
technique of CH 4 differential imaging. 
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Table 5 

Polynomial coefficients for plotted fits 



P(x 


) 




RMS 




ci 




c 3 


C4 














M K 


a (2M) 


SpT 


0.38 


6.27861c+l 


-1.47407c+l 


1.54509c+0 


-7.424180-2 


1.63124c-3 


-1.25074c- 


-5 










M K 


(U) 


SpT 


0.40 


8.14626c-l 


2.95440e+0 


-3.89554e-l 


2.68071o-2 


-8.868850-4 


1.14139e- 


-5 










Mj 


(2M) 


SpT 


0.36 


8.94012e+2 


-3.98105e+2 


7.57222e+l 


-7.86911e+0 


4.82120e-l 


-1.73848c- 


-2 


3.41146c- 


-4 


-2.80824c- 


■6 


Mj 


(U) 


SpT 


0.30 


5.04642o+l 


-3.13411e+l 


9.06701o+0 


-1.30526o+0 


1.03572e-l 


-4.58399c- 


-3 


1.05811e- 


-4 


-9.91110c- 


■7 


M Jc 


SpT 


0.37 


7.22089e+l 


-1.58296e+l 


1.56038e+0 


-6.49719e-2 


1.04398e-3 


-2.49821e- 


-6 










M z 


(U) 


SpT 


0.29 


4.99447o+l 


-3.08010e+l 


8.96822e+0 


-1.29357e+0 


1.02898e-l 


-4.57019e- 


-3 


1.05950e- 


-4 


-9.97226c- 


■7 


M Ic 


Ic-J 


0.67 


1.17458e+3 


-1.22891e+3 


5.00292o+2 


-9.70242e+l 


8.86072e+0 


-2.97002e- 


-1 










Mj 


(2M) 


Ic-J 


0.63 


1.52199c+3 


-1.69336c+3 


7.41385o+2 


-1.58261c+2 


1.64657c+l 


-6.66978c- 


-1 











a SpT — i for Mi, — j + 10 for Lj and — z + 19 for Tz spectral types on the Kirkpatrick ct al. (1999) system for M and L dwarfs, and the Burgasser et al. 
(2002b) system for T dwarfs. 

6 P{x) — CO + CI X + C2 X 2 ... 
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